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Abstract 

In this paper we investigate a stochastic model for an economic 
game. To describe this model we have used a Wiener process, as the 
noise has a stabilization effect. The dynamics are studied in terms 
of stochastic stability in the stationary state, by constructing the 
Lyapunov exponent, depending on the parameters that describe the 
model. The numerical simulation that we did justifies the theoretical 
results. 
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1 Introduction. 

Stochastic modeling plays an important role in many branches of sci- 
ence. In many practical situations, perturbations are expressed in terms of 
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white noise, modeled by brownian motion. The behavior of a deterministic 
dynamical system which is disturbed by noise may be modeled by a stochastic 
differential equation (SDE), [7]. The stochastic stability has been introduced 
by Bertram and Sarachik and is characterized by the negativeness of Lya- 
punov exponents. In general, it is not possible to determine this exponents 
explicitly. Many numerical approaches have been proposed, which generally 
used the simulation of the stochastic trajectories [3]. In the present paper, we 
study a stochastic dynamical system that are used in economy, in describing 
a Counot duopoly game. 

In 1838, Cournot introduced the first formal theory of oligopoly, which 
treated the case of naive expectations, where each player assumes the last 
values taken by the competitors without estimation of their future reactions 
[4]. Recently, a lot of articles have shown that the Cournot model may lead 
to a cyclic or chaotic behavior [2], [3], [8], [10], [11], [12]. Also, in Rosser 
reviews the development of the theory of complex oligopoly dynamics. 

In the present paper we have studied a stochastic Cournot economic game. 
In Section 2 we present the Lyapunov exponent and stability in stochastic 
2d dynamical structures. Section 3 studies the Lyapunov exponent for an 
economic game with stochastic dynamics. Some numerical simulations are 
given in Section 4. Finally, Section 5 draws some conclusions. 

2 The Lyapunov exponent and stability in 
stochastic 2d dynamical structures. 

Let {Q, 3", y) be a probability space [7|. It is assumed that the a— algebra 
3" is a filtration that is, 3" is generated by a family of cr— algebra 3^t{t > 0) 
such that 

3'sC3'tC3', Vs<t,s,tG/, 

where I =[0,T],T e (0,oo). 

Let {x(t) = {xi(t),X2(t))}t>o be a stochastic process. The system of Ito 
equations: 

dxi{t, u) = fi{t, x{t, uj))dt + gi{x{t, uj))dw{t, uj),i = 1,2, (1) 
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with the initial condition x(0) = xq is written as: 



Xi{t,io) = Xio{io) + fi{x{s,uj))ds + gi{x{s,uj))dw{s,u),i = 1,2, (2) 
Jo Jo 

for almost all u; G ^2 and for each t > 0, where fi{x) are drift functions, 
gi{x) are diffusion functions, fi{x{s))ds, i = 1,2 are Riemann integrals and 
/o 9i{-'-{s))dw{s) are Ito integrals. It is assumed that and gi, i = 1,2 satisfy 
the conditions of existence of solution for this SDE with initial condition 
x(0) = Co G IR". 

Let Xq = {xio, X20) € IR^ be a solution of the system: 

Mxo)=0,t = l,2. (3) 

The functions gi,i = 1,2 are chosen so that: 

gi{xo) = 0,z = 1,2. 

In what follows, we consider: 

2 

9i{^) = bij{xj — Xoj), i = 1,2, 

where bij G = 1,2. 

The linearized system of ([2]) in xq, is given by: 



where 



X{t) = [ AX{s)ds + [ BX{s)dw{s) 
Jo Jo 



Xit)=[ V^=f'''' ^^0,5-^^" 



y{t,uj) y ' \ 021 a22 / ' V ^21 &22 



The Oseledec multiplicative ergodic theorem [9] asserts the existence of 2 
non-random Lyapunov exponents A2 < Ai = A. The top Lyapunov exponent 
is given by: 

A = lim sup log \J x{tY + y{tY. 



3 



Applying the change to polar coordinates: 

x{t) = r{t)cose{t),y{t) = r{t)sine{t) 
by writing the ltd formula for 

1 y 

hi{x,y) = - log{x^ + y^) = \og{r), h2{x,y) = arctg{-) = 9. 

2 X 

we get: 

Proposition 1 /5/. The formulas 
log (^) = f\Ms))+kqMs)r-q2i9{s)r)ds+ f qMs))dw{s), (4) 



eit) = 9i0)+ I q,{eis))-q2{e{s)qMs))ds+ I qMs))dw{s), (5) 
hold, where 



qi{0) = auCos'^{6) + (ai2 + a2i)cosd sinO + a22sin'^0, 
q^{B) = bucos'^{0) + {bi2 + b2i)cos6sm9 + b22sin'^0, 
q3{9) = a2iCos'^{9) + (022 — aii)cos6 smO — ai2sin^9, 
qA^O) = b2iCos^{9) + (^22 — bii)cos9sm.9 — bi2sin^9. 

As the expectation of the Ito stochastic integral is null 

E [ q2{9{s))dw{s) = 0, 
Jo 

the Lyapunov exponent is given by: 

A= hm ilog (^) = hm ^-E [\qMs)) + kqMsW -q2{0{s)))ds. 
Applying the Oseledec theorem, if r(t) is ergodic, we get: 

A = l\qi{9) + ^(g3(^)' - q2{9)))p{9)d9, 
where p{9) is the probability distribution of the process 9. 



(6) 
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An approximation of this distribution is calculated by solving the Fokker- 
Planck equation. 

The Fokker-Planck (FPE) equation associated with equation ([5]) for p = 
pit, 9) is 

f + - q2{e)q4{e))p) - l^2i^m'p) = 0. (7) 

From ([7]), it results that the solution p{6) of the FPE is solution of the 
following first order equation: 

{-qm+qiio)qm+q2{e)g,{e))p{e)+^q,{e)y{e)=po, (8) 

dp 

where p'{9) = — and 
d9 

q^ie) = -{bi2 + b2i)sin2e - (622 - bn)cos2e. 

Proposition 2 f^. If qi{9) 7^ 0, the solution of the equation ^ is given 
by: 

where k is determined by the normality condition 

/ p{e)d9 = 1 

Jo 

and 

D{2n) - 1 
lo^ D{u)du 

The function D is given by: 

Die) = exp(-2 r g3H-g2Hg4H -g4Hg5H^^^ 
Jo qiiur 

A numerical solution of the phase distribution could be performed by a 
simple backward difference scheme. 
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We consider N e IR+, ^ = 

qi{i) = ail cos^(i/i) + (ai2 + 021) cos{ih) sm{ih) + 022 sin^(i/i), 
g2(^) = ^11 cos^(i/i) + (&12 + ^21) cos{ih) sm{ih) + 622 sin^(i/i), 
g3(i) = a2i cos^(i/i) + (022 — Oil) cos{ih) sm{ih) — ai2sin'^{ih), 
q4{i) — 621 cos^(i/i) + (622 — &11) cos{ih) sm{ih) — bu sin^(i/i), 
Qbii) = -(^'12 + hi) sin(2i/i) - (622 - bu)cos{2ih),i = 0, N 

The function p{i), i = 0, N is given by tlic following relations: 
where 



The Lyapunov exponent is A = A(A^), where 

^ 1 

A(iV) = + -(g4(^f - ?2m)pW/.. 

Proposition 3 // the matrix B is given by: 

bu = a, bi2 = -/3, 621 = ^22 = a 

A; 1 1 1 

p(0) = ^exp{^((a2i - ai2 - cxP)e + -{an - 022) cos 26*+ -(021- 

— 012) sin 2^)} 



k- 



/Q^''exp{-^((a2i-ai2-Q;/3)6'+^(aii-a22) cos 26'+^ (021-012) sin26')d6' 
A = ^(aii + 022 + - a^) + ^(aii - 022)02 + ^(021 + 012)^2, 

where 

C2= / cos{2e)p{e)d0, S2= / sin{29)p{0)de. 
Jo Jo 



3 The Lyapunov exponent for an economic 
game with stochastic dynamics. 

Two firms enter the market with a homogenous consumption product. 
The elements which describe the model are: the quantities which enter the 
market from the two firms Xi > 0, i = 1,2; the inverse demand function 
p : — s> ]R+ ( p is a derivable function with p' (x) < 0, limp{x) = 0, 

x~*ai 

limp{x) = bi, (oi G M, 6i G M); the cost functions Cj : M+ M+ ( Cj are 
derivable functions with Q (xj) > 0, > 0, i = 1, 2 ). 

In our study we consider p{x) = —,x > and Ci{xi) = CiXi + di,i = 1, 2. 

The mathematical model of the stochastic dynamic economic game is 
described by the stochastic system of equations: 

X2(s) /"* 

xi{t) = xi{0)+ki {- — — rT^-Ci)ds+ {bnXi{s) + bi2X2is) + -fj^)dw{s) 

Jo [Xl{S)+X2[S)) Jq 

X2{t) = X0)+k2 I — -—T-C2)ds+ {b2iXi{s) + b22X2{s)+'j2)dw{s) 

Jo [Xl{s)+X2{s)y Jo 

(9) 

where bij G IR, i,j = 1,2, ki > 0,k2> 0, Xi{t) = Xi(t,u!), i = 1,2. 

_ bnC2 + &12C1 _ b2iC2 + &22C1 

For bij = 0, i, j = 1, 2 model Qj is reduced to the classical model of the 
economic game [2], [H]- 

The system of stochastic equations ([9D, (SDE), has the form ([2]) from 
section 2, where: 

fi{xi,X2) = 7 — 7^^— - Ci,gi{xi,X2) = bnXi + bi2a;2 + 7i, 

[Xi + X2Y 

f2{Xl,X2) = 7 ■ rij - C2,g2{Xl,X2) = b2lXi + b22X2 + 72- 

[Xi + X2Y 

Applying the results from section 2, we have: 
Proposition 4 (i) The stationary state of (SDE) ^ is given by: 

C2 Ci 



^10 — 7 ; r?) ^20 



(Cl+C2)2'-^^'^ (Ci+C2)2' 
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(a) The elements of the matrix A, which characterize linearized equation 
^ in (xio, 2:20) are: 

au = -2fciCi(ci + C2), ai2 = -ki{cl - cj) 
a2i = k2{cl - cl),a22 = -2k2C2{ci + C2); 

(Hi) The roots of the characteristic equation: 

- («11 + «22)/W + «11«22 - «12«21 = (10) 

have the real part: 

^e(/ii_2) = -{kici + k2C2){ci + C2); 

(iv) If bii = a, bi2 = —[3, 621 = [3, 622 = a, /? 7^ 0, then the Lyapunov 
coefficient of (SDE) ^ is: 

X = -{kiCi + k2C2){ci + C2) + - a^) - {kiCi - k2C2){ci + 02)02+ 

1 ^ (11) 

+ -{k2-k^){cl-cl)E2 

where 

r.27r r2iT 



D2= cos{2d)p{e)de,E2= / sin{2e)p{d)de 
Jo Jo 



and 



p{e) = kg{e),k 



lo 

2exp{-2, 



g{9) = ^exp{^((A;i + k2){cl - cl) + a(3)9 - {kiCi - k2C2){ci + 02) cos(2^) + 
1 



+ -(k, + k2)icl-cl)smi2e)}. 



4 Numerical Simulations. 



We have done the numerical simulations using a program in Maple 12. 
For Cl = 0.2, C2 = 2, ki = 0.2, ^2 = 0.4, /3 = 2, in figure 1 is displayed 
(a, A(a)), where A(a) is given by ffTTl) . For a G (—00, —1.2) U (1.1, cxd), the 
Lyapunov exponent is negative, then (SDE) has an asymptotically stable 
stationary state. For a G (—1.2, 1.1), the Lyapunov exponent is positive and 
(SDE) has an asymptotically unstable stationary state. 
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Fig 1. (a,A(a)) 




If is a real parameter and a = 2, the figure 2 shows the behavior of the 
top Lyapunov exponent as a function of (3: {(3, A(/?)). 



Fig 2. {p,\m 




For (3 G (— oo, — 2.6) U (2.6, oo) the Lyapunov exponent is positive and 
(SDK) has an asymptotically unstable stationary state. For /?(— 2.6, 2.6) 
the Lyapunov exponent is negative and (SDE) has an asymptotically stable 
stationary state. 
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The Euler second order scheme for (SDE) (2) is given by: 

X2(n) 



Xi{n + 1) = Xi(n) + h 



+ X2{n)) 



2 ci )+ (611X1 (n) + bi2a;2(ri) + 7i) 



■ G{n) + f>ii(feiixi(n) + h2X2in) + 7i)^^^^^7^ + 



2xi{n)x2{n) 
{xi{n) + X2{n))^ 



X2{n) \ n / \ 7 / \ N Xi(n)x2(n) , /i^ 

- ci + (feiiXiH + 6i2X2(n) + 7i) . . ^ \ \.J — + 
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{xi{n) + X2(n))2 
2z2(n) 
(xi(n) + X2{n)y 



(xi(n) + X2(n))3 2 



)(6iiXi(n) + 6i2X2(n) + 71) 



hGjn) 



X2{n + 1) = X2{n) + h 



Xi{n) 



{xi{n) + X2{n)y 

G{nf - h 



■ G{n) + 622(621X1(72) + h22X2{n) + 73) 



C2 ) + (6212:1 (n) + 622X2(71) + 72) 

2x1 (72)^2 (n) 
^ {xi{n) + X2{n)Y 



xAn) \ /, / X , /X X Xi(7l)x2(72) , /i^ 

^ ^ - C2 + (621X1(72) + 622X2(77) + 72) . / x \\xJ — + 



(621- 



(^1(77) + ^2(77))^ 
2X1(77) 



(^1(77) + X2(77))3 2 



(Xi(77) + X2{n)Y 



) (621X1 (77) + 622X2(77) + 72) 



/7G(n) 



where G'(72) = 77;((77 + l)h) — w{nh), n = 1,2,..., and Xj(77) = Xi{nh,uj), 
7 = 1,2. 

In figures 3 and 4 are displayed the orbits: (77, xi(77, cj)) for (SDE) and 
(77,xi(77)) for (ODE): 



Fig 3. (77, xi(77, uj)) 



Fig 4. (77, Xi(77)) 
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In figures 5 and 6 are displayed the orbits: {n,X2{n,uj)) for (SDE) and 
(n,X2(n)) for (ODE): 



Fig 5. {n,X2{n,uj)) 



Fig 6. (n, X2{n)) 



20O \ 400 600 KOO 1000 



200 400 600 800 1000 



In figures 7 and 8 are displayed the orbits: {xi{n,u!),X2{n,uj)) for (SDE) 
and {xi{n),X2{n)) for (ODE): 



Fig 7. {xi{n,u),X2{n,uj)) Fig 8. {xi{n) , X2{n)) 



0.36 0.37 0.35 0.3? 0.40 :0.41 0.42 0.43 0.44-0.45 



0-400 0,402 0,404 0,406 0,40S 0,410 0,412 



5 Conclusions. 

In the present paper we investigate an economic game with stochastic 
dynamics. We focus on a particular game and determine the Lyapunov expo- 
nent for the stochastic system of equations that describes the mathematical 
model. The calculation of the top Lyapunov exponent enables us to decide 
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whether a stochastic system is stable or not. Using a program in Maple 12, 
we display the Lyapunov exponent and the system orbits. 

References 

[1] Arnold, L., Random dynamical systems, Springer Monographs in Math- 
ematics, Springer- Verlag, Berlin, 1998. 

[2] Bundau O., Neamtu M., Opris D., Rent seeking games with tax evasion, 
Pannonian Applied Mathematical Meetings, B.A.M-CX, 2007; 2295- 
2315: 121-128; arXivrmath.DS. 0706.0664vl, 6 June, 2007. 

[3] Bischi G.I., Naimzada A., Global analysis of a dynamic duopoly game 
with bounded rationality. In Advances in Dynamic Games and Applica- 
tions, Vol. 5., Birkhaur, Boston, 1999. 

[4] Cournot A., Researches into the principles of the theory of wealth. Engl, 
trans.. Chapter VII, Irwin Paper Back Classics in Economics, 1963. 

[5] Jedrzejewski, P., Brochard, D., Lyapunov exponents and stability in 
stochastic dynamical structures. 

[6] Hu. B.Y., Mahommcd, S.E., Yan, P., Discrete-time approximation of 
stochastic delay equations. The Annals of Probability, vol. 32, Nr lA 
(2004), 265-314. 

[7] Kloeden, P.E., Platen, E., Numerical solution of stochastic differential 
equations. Springer- Verlag, 1992. 

[8] Mircea G., Neamtu M., Opri§ D.: Dynamical systems from economy, 
mechanics, biology described by differential equations with time delay 
(in Romanian), Mirton Publishing house, (2003). 

[9] Oseledec, V.I., A multiplicative Ergodic theorem, Lyapunov character- 
istic numbers for dynamical systems, Trans. Moscow Math. Soc. 1968, 
no.l9, 197-231. 

[10] Puu T., Chaos in duopoly pricing. Chaos, Solitons and Fractals 1991; 1: 
573-581. 



12 



[11] Puu T., The Chaotic duopolists revisited. J. of Economic Behavior and 
Organization Vol. 1998; 37: 385-394. 

[12] Puu T., Attractors, bifurcations and chaos: Nonlinear phenomena in 
economics, Springer- Verlag, 2000. 

[13] Rosser J.B., The development of complex oligopoly dynamics theory. In 
Text Book Oligopoly Dynamics: Models and Tools Springer -Verleg, 
2002. 

[14] Schurz, H., Moment contractivity and stability exponents of nonlinear 
stochastic dynamical systems, IMA Print Series [^1656, 1999. 



13 



